High-density QCD and Cosmic Ray Air Showers 
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We discuss particle production in the high-energy, small-a; limit of QCD where the gluon density 
of hadrons is expected to become nonperturbatively large. Strong modifications of the phase-space 
distribution of produced particles as compared to leading-twist models are predicted which reflect in 
the properties of cosmic ray induced air showers in the atmosphere. Assuming hadronic primaries, 
our results suggest a light composition near GZK cutoff energies. We also show that cosmic ray data 
discriminate among various QCD evolution scenarios for the rate of increase of the gluon density at 
small x, such as fixed-coupling and running-coupling BFKL evolution. There are clear indications for 
a slower growth of the gluon density as compared to RHIC and HERA, due e.g. to running-coupling 
effects. 
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Today, the cosmic ray energy spectrum has been measured up to energies near the GZK cutoff, E w 10 11 GeV 0,0. 
These energies by far exceed those reached by terrestrial accelerators. Thus, air showers induced in our atmosphere 
"^s 1 present a unique opportunity to probe high-energy QCD at very small light-cone momentum fractions x, i.e. in the 
regime of nonperturbatively large gluon densities. The physics of gluon saturation is therefore expected to play a 
significant role for the properties of extensive air showers, and thus for the determination of the nature of the highest 
energy cosmic rays. We refer to refs. 0, 0| for discussions regarding the relevance of high energy QCD interactions 
^ . for air showers and composition analysis. 

The high-energy limit of hadron scattering from a nucleus can be adressed from two complementary views. In the 
frame where the nucleus is at rest the partons up to the "black-body" resolution scale pt(s) interact with the target 
00 ■ with (nearly) the geometric cross section of 2-kR 2 a . Hence, in this limit the projectile wave function is resolved at 
a virtuality of ~ p\ which is much larger than any soft scale such as Aqcd- In this frame, the process of leading 
hadron production corresponds to releasing the resolved partons from the projectile wave function. The partons then 
fragment with large transverse momenta ~ pt and essentially independently since their coherence was completely lost 
in the propagation through the black body. In the case of 7* A scattering one is able to make nearly model independent 
O |' predictions for the leading hadron spectrum [j| which differ drastically from the DGLAP leading-twist limit. 
' , On the other hand, one could discuss the high density limit in the infinite momentum frame. Indeed, the wave 
function of a fast hadron (or nucleus) exhibits a large number of gluons at small x. The density of gluons is expected to 
saturate when it becomes of order l/a s Q . The density of gluons per unit of transverse area at saturation is denoted 
by Q 2 , the so-called saturation momentum. This provides an intrinsic momentum scale @ which grows with atomic 
number (for nuclei) and with rapidity, due to continued gluon radiation as phase space grows. For sufficiently high 
energies and/or large nuclei, Q s can become much larger than Aqcd and so weak coupling methods are applicable. 
Nevertheless, the well-known leading-twist pQCD can not be used precisely because of the fact that the density of 
gluons is large; rather, scattering amplitudes have to be resummed to all orders in the density. When probed at a scale 
below Q s , cross sections approach their geometrical limit over a large range of impact parameters, while far above Q s 
one deals with the dilute regime where they can be approximated by the known leading-twist pQCD expressions. 

The target nucleus, when seen from the projectile fragmentation region, is characterized by a large saturation 
momentum. Its precise value can not be computed from first principles at present but model studies of deep inelastic 
scattering (DIS) at HERA suggest 

Ql{x) ~ Ql (x /x) X (1) 

with xq a reference point and an intercept A » 0.3 @. The initial condition at xo accounts for the growth of Q s with 
the number of valence quarks; for example, near the rest frame of the nucleus one might fix Qq °c A 1 ! 3 log A, with a 
proportionality constant of order Aqcd Q- [We remark that for realistic nuclei Q s does of course also depend on the 
impact parameter, which we don't spell out explicitly] 

The above scaling relation can be obtained from the fixed coupling BFKL evolution equation for the scattering 
amplitude of a small dipole. The BFKL equation is a linear QCD evolution equation which can not be applied in 
the high-density regime. Nevertheless, one can evolve the wave function of the target in rapidity y = log 11 x and 
ask when the dipole scattering amplitude becomes of order one, which leads to Q 2 (y) = Q 2 (yo) exp ca s y jjjl, with 
a s = a s N c /n and c « 4.84 a constant. Hence, LO fixed-coupling BFKL evolution predicts A' = ca s of order one, a 
few times larger than the fit (QJ to HERA phenomenology. A NLO BFKL analysis corrects this discrepancy and leads 



2 



to A' much closer to the phenomenological value A similar observation is made in ^lj where both log(l/x) and 
logQ effects were considered. The approach from the leading-twist regime to the black-body limit (BBL) has also 
been studied using the DGLAP evolution equation in ref. |l2l |. 

On the other hand one could also consider BFKL evolution with ad-hoc one- loop running of the coupling [flL IT3|: 

^QCd); 



a s {Q 2 s ) = & /log(Qg(a;)/A| ir! r»'), which leads to 



Ql(y) = a qcd exp y/2b c{y- 
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with 2&ocyo = l°g 2 (Qo/AQCD)- Insisting that be valid at least in the y — ► limit again provides us with a 
phenomenological value for the constant c in terms of the saturation momentum at y = 0. The form leads to a 
notably slower growth of Q s at high energy. Specifically, for central proton-nitrogen collisions at RHIC, LHC and 
GZK-cutoff energies (total rapidity y — 10.7, 17.3 and 26.0) the saturation momentum of the nucleus in the rest frame 
of the projectile hadron is Q s = 1.4, 4.5, 19.2 GeV for fixed coupling evolution, while for running coupling evolution it 
is Q s = 1.1, 2.4, 5.9 GeV, respectively. Clearly, cosmic ray interactions in our atmosphere offer a realistic opportunity 
for observing this effect. 

The dominant process for fast particle production (xf > 0.1) is scattering of quarks from the incident dilute 
projectile on the dense target. At high energies (i.e. in the eikonal approximation where p + is conserved), the 
transverse momentum distribution of quarks is given by the correlation function of two Wilson lines V running along 
the light cone at transverse separation r t (in the amplitude and its complex conjugate), 
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Here, the convention is that the incident hadron has positive rapidity, i.e. the large component of its light-cone 
momentum is P + , and that of the incoming quark is p + = xP + (q + for the outgoing quark) . The two-point function 
has to be evaluated in the background field of the target nucleus. A relatively simple closed expression can be 
obtained ^4[ m the McLerran-Venugopalan model of the small- a; gluon distribution of the dense target Q . In that 
model, the small- a; gluons are described as a stochastic classical non-abelian Yang-Mills field which is averaged over 
with a Gaussian distribution. The qA cross section is then given by 0] 
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This expression is valid to leading order in a s (tree level), but to all orders in Q s since it resums any number of 
scatterings of the quark in the strong field of the nucleus. The saturation momentum Q s , as introduced in eq. (@J, is 
related to x, the total color charge density squared (per unit area) from the nucleus integrated up to the rapidity y 
of the probe (i.e. the projectile quark), by Ql — 47r 2 a^ \ ~ l)/N c . In the low-density limit, x ls related to the 
ordinary leading-twist gluon distribution function of the nucleus, see for example [l^ . 

The integrals over pt in eq. (QJ are cut off in the infrared by some cutoff A, which we assume is of order Aqcd- 
At large transverse momentum, the first exponential in <QJ can be expanded order by order to generate the usual 
power series in l/qf with the leading term corresponding to the perturbative one-gluon i-channel exchange 

contribution to oo — » qg scattering. On the other hand, for Q s ~ q t 3> A one obtains in the leading logarithmic 
approximation 
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Qs, and hence the transverse 



This approximation reproduces the behavior of the full expression iQJ about 
momentum integrated cross section reasonably well. It is useful when the cutoff A <C Q s , that is, when color 
neutrality is enforced on distance scales of order 1/A> 1/Q S - If, however, color neutrality in the target nucleus 
occurs over distances of order 1/Q B E3l then A ~ Q s and one would have to go beyond the leading- logarithmic 
approximation; work along those lines is in progress. 

Consider the probability of inelastic scattering (i.e. with color exchange) to small transverse momentum. This is 
given by expression 1)415(1 . integrated from q t = to q t = A: 



d 2 q t C(q t ) 
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Q 2 log Q s /A 



(6) 
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Here, we have written only the leading term in A 2 /Q 2 , neglecting subleading power-corrections and exponentially 
suppressed contributions. Hence, soft forward inelastic scattering is power-suppressed in the black-body limit. This 
steepens the longitudinal distribution dN/dxF of leading particles since partons with large relative momenta fragment 
independently (3.Hfl|. 

We now turn to gluon bremsstrahlung which dominates particle production at xf < 0.1. Gluon radiation with 
transverse momentum q t ~ Q s in high-energy hadron-nucleus collisions has been discussed in detail in [111 El; nere 
we employ the latter ansatz for the fusion of gluon ladders: 

2 

it 

E ir= 47T ^^^2 f dk 2 t a s {k 2 )<f> h {x u k 2 )<}> A {x 2 ,{q t -k t ) 2 ) , (7) 
dPq N£ - 1 qf J 

where 4>(x,Q 2 ) denotes the unintegrated gluon distribution function of the projectile hadron or target nucleus, re- 
spectively. It is related to the gluon density by 

Q 2 

xg(x,Q 2 ) = J dk 2 4>{x,k 2 ) . (8) 

The small- x gluon density has been investigated recently within RG-improved QCD evolution equations by Ciafaloni 
et al. [Tl] |. It predicts the onset of nonlinear (saturation) effects at transverse momenta similar to our values for Q s {x). 
Their gluon densities can be incorporated into phenomenological applications in the future. 

For practical reasons, we presently employ the simpler Kharzeev-Levin ansatz for the infrared-finite gluon densities 

x g(x, Q 2 ) oc — min(Q 2 , Q 2 s {x)) (1 - xf , (9) 

with a s evaluated at max(Q 2 , Q 2 ). In our approach, the absolute normalization is determined by energy conservation. 
The number of produced gluons behaves as dN / dq 2 ~ l/g^ at large transverse momentum but flattens to ~ 1/qf at 
qt < max(Qf(x),Q^(x)) and finally approaches a constant in the qt — ► limit. Despite the rather qualitative nature 
of this ansatz, the main feature is that saturation effects provide an intrinsic semi-hard scale Q s for partonic processes, 
thus eliminating the infrared divergences of leading-twist perturbation theory and, at the same time, the need for 
matching to some purely phenomenological models for low-g* particle production. In fact, the effects discussed here 
do not depend on details of our approach but follow generically from the fact that the average transverse momentum 
(~ Qs) gets large as the BBL at high energy is approached. 

For the application to cosmic-ray induced airshowers we have developed a Monte-Carlo algorithm which generates 
complete parton configurations as described above. The partons are then connected by strings. This accounts for the 
"absorption" of almost collinear gluons on the string, while independent fragmentation is reproduced when relative 
transverse momenta of partons are large. The fragmentation is performed via the Lund scheme as implemented in 
PYTHIA Wk. T his model is linked to a standard pQCD event generator commonly used in air-shower computations 
(Sibyll v2.1 |2flj) which handles low-energy and peripheral collisions where the saturation momentum of the nucleus 
is not sufficiently large. Finally, the hadron-nucleus collision models are embedded into the cascade equations which 
solve for the longitudinal profile of the airshower [24] . Details of our Monte-Carlo implementation will be published 
elsewhere. 

Fluorescence detectors measure the number of charged particles (mostly e ) at a given atmospheric depth X which 
is given by the integral of the atmospheric density along the shower axis, X = J ds p(s). The position of the maximum 
defines X max which increases monotonically with the energy of the primary. Note that for nuclei the primary energy 
is shared by all of its nucleons and so A max also depends on the mass number: at fixed E, heavier primaries lead to 
smaller "penetration depth" X max , 

In FigureQwe compare the predictions of the leading-twist pQCD model Sibyll for proton and iron induced showers 
to the saturation model (BBL, for proton primaries only) with running and fixed coupling BFKL evolution of Q s , 
respectively, and to preliminary Hires stereo data p|. In the saturation limit, showers do not penetrate as deeply 
into the atmosphere. This is due to the "break-up" of the projectile's coherence ^| together with the suppression of 
forward parton scattering (for central collisions). The comparison to the data suggests a light composition at those 
energies. 

Also, contrary to present accelerator experiments, a clear difference between running-coupling and fixed-coupling 
BFKL evolution of the saturation momentum is apparent in this observable. The discrepancy between those evolution 
scenarios at the highest energies is strongly amplified by subsequent hadronic collisions in the cosmic ray cascade since 
it determines the fraction of events that occur close to the black-body limit (averaged over all impact parameters). 
Thus, assuming hadronic primaries, the extremely rapid growth of Q s obtained for fixed coupling evolution is at 
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Figure 1: X max as a function of primary energy. 

variance with the Hires data, as it would require hadrons lighter than protons. This is due to a too strong suppression 
of leading hadron production over a large range of impact parameters at high energies. At lower energies, of course, 
the two evolution scenarios predict similar saturation scales and so can not be distinguished as reliably by present 
collider experiments. 

Finally, we remark that our results for running-coupling evolution coincide with those of another popular hadronic 
model, QGSJET f^]. Due to the absence of an ad-hoc qt cutoff for pQCD interactions in QGSJET, that model 
needs to assume a rather flat gluon density ~ I/a/ 2 ' a t small x in order not to overestimate multiplicities at collider 
energies fj|. In our approach, on the other hand, the increase of the multiplicity and of the typical transverse momenta 
with energy is controlled by the saturation mechanism and the corresponding evolution of the gluon density. 

In conclusion, we have shown that at energies near the GZK cutoff QCD evolution scenarios differ drastically in 
their predictions for the scale Q s where gluon densities become nonperturbatively large (the BBL). Assuming hadronic 
primaries, properties of induced air showers suggest a moderately rapid growth of the gluon density with decreasing 
x (over a wide range of x), corresponding to a more slowly increasing black-body scale Q s with log(l/x). At least 
qualitatively, this agrees with a number of recent studies u sing extensions of DGLAP, resummations of log x and log Q 
effects, and running-coupling BFKL evolution |ir]L fllL fl^J . fl7^ | . 

Fixed-coupling BFKL evolution, which is able to fit DIS at HERA Q and deuteron-gold multiplicities at RHIC [2lj . 
predicts extremely large values for Q s at GZK energies. From our analysis we conclude that its applicability is limited 
to much lower energies. We stress that measurements of hadron spectra in pA collisions at the LHC over a broad 
region of phase space (in particular, in the forward region) could majorly advance our knowledge of high-density QCD 
and, in turn, help us understand the nature, composition and perhaps the origin of the highest-energy cosmic rays. 
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